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EXPERIMENTAL 

Precision  measurements  of  volume  versus  pressure  for  one  meter  long 
potassium  specimens  have  been  successfully  carried  out  up  to  a  hydrostatic 
pressure  of  7  kbar.  The  length  change  of  the  specimen  has  been  measured  to 

o 

an  accuracy  of  better  than  500  A  using  a  He-Ne  laser  interferometer  in  a 
temperature  controlled  environment  with  temperature  variations  ±0.001°C  at 
28.58°C. 

Extensive  work  was  required  to  grow  these  large  single  crystals  of  a 
reactive  material  and  considerable  work  was  needed  to  develop  techniques  for 
handling  these  large,  reactive,  compliant  and  weak  crystals.  The  successful 
growth  of  this  crystal  can  only  be  described  as  a  tour  de  force.  It  took  us 
much  longer  than  expected  because  it  was  much  harder  than  expected. 

Potassium  immersed  in  a  mineral  oil  environment  inside  a  5  gallon  pail  was 
purchased  from  Mine  Safety  Appliance  Corporation,  Evans  City,  PA,  with  a  99.95% 
purity.  It  was  transferred  from  the  5  gallon  pail  into  an  extrusion  vessel 
specially  made  for  handling  potassium.  Since  potassium,  like  other  alkali 
metals,  is  very  soft  and  oxidizes  rapidly  in  moisture  and  air,  transfer  was 
made  directly  into  the  extrusion  vessel  filled  with  liquid  hexane  which  had 
been  previously  purified  and  dehydrated  by  keeping  it  in  the  container  with  a 
large  amount  of  potassium  chips  under  slight  pressure  of  ultra  high  purity 
argon  gas.  Then  the  extrusion  was  carried  out  into  a  pyrex  glass  tube  which 
was  also  filled  with  pre-purified  liquid  hexane  to  a  desired  length  (about 
130  cm  long)  using  a  7.9  mm  die.  This  extruded  potassium  was  zone-refined 
several  times  to  enhance  purity  and  anneal -out  defects  and  some  of  the  dis¬ 
locations.  It  was  then  finally  transferred  to  our  length  measurement  pressure 
vessel.  Liquid  hexane,  purified  as  described  above,  was  used  as  a  pressurizing 
medium. 
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Our  data  have  been  analyzed  by  a  special  non-linear  least  square  fit  of 
various  equations  of  state  to  obtain  the  isothermal  bulk  modulus  and  its 
pressure  derivatives  at  atmospheric  pressure. 

The  first  order  Murnaghan  equation  (ME-j)  which  assumes  that  all  pressure 
derivatives  of  the  bulk  modulus  at  zero  pressure  beyond  the  first  are  zero,  the 
second  order  Murnaghan  equation  (ME2),  two  parameter  Birch  equation  (BE-j ) , 
three  parameter  Birch  equation  (BE2),  and  Keene  equation  yield  the  follow¬ 
ing  results  listed  in  Table  I  below. 


Table  I. 

Bulk  modulus  and 
at  28.58°C. 

its  pressure  derivatives 

for  potassium 

Equation 

Used 

B0 

(kbar) 

Bo' 

B0"  1 

(kbar)  1 

MEi 

31.01±0.02 

3.905+0.009 

me2 

30. 87+0.01 

4.094+0.018 

-0.0747t0.0068 

BE1 

30.85±0.01 

4.144+0.003 

-0.131t0.002 

CVJ 

LU 

CD 

30.85±0.02 

4.151 ±0 . 021 

-0.135t0.01 

KE 

30.84  0.02 

4.159+0.024 

-0.145i0.014 

In  the  above  table,  B  ,  B  '  and  B  "  are  the  values  of  bulk  modulus,  its 

0  0  0 

first  pressure  derivative,  and  its  second  pressure  derivative,  taken  at  atmos¬ 
pheric  pressure,  respectively.  The  above  data  are  to  be  compared  with 
Swenson's  values  [C.  E.  Monfort  and  C.  A.  Swenson,  J.  Phys.  Chem.  Solids  26, 
291  (1965)]  obtained  by  the  piston-cylinder  method  and  Smith's  ultrasonic 
measurements  [P.  A.  Smith  and  C.  B.  Smith,  J.  Phys.  Chem.  Solids  26,  279 
(1965)].  %heir  data  were  fitted  by  ME1  which  is  based  on  the  assumption 
that  B  =  BQ  +  B0'P,  and  yield 

B(kbars)  =  30.8  +  3.85  P  (Swenson) 

B(kbars)  =  30.9  +  3.98  P  (Smith) 
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The  manner  of  fitting  is  itself  very  important;  note  how  B0'  of  our  data  changes 
when  ME,  is  used  rather  than  ME,.  Thus  Swenson's  B  '  would  be  increased  from 
3.85  if  MEg  or  other  higher  order  equations  were  used.  It  is  quite  exciting 
to  see  that  our  values  for  BQ  are  in  excellent  agreement  with  those  obtained 
by  Swenson  and  Smith  and  that  the  BQ'  values  are  close.  It  should  be  noted 
that  various  theories  of  binding  show  that  the  Maclauren  series  for  B(P)  in¬ 
volves  alternating  signs  for  the  derivatives,  with  BQ'  positive,  B0"  negative, 
etc.  Hence  the  use  of  ME-j  will  give  a  lower  bound  on  BQ'  and  an  upper  bound 
on  Bq,  because  of  the  neglect  of  the  negative  B0"P2/2  term.  Similarly,  the 
use  of  ME2  should  give  a  lower  bound  on  | B0" ) ,  an  upper  bound  on  BQ'  and  a 
lower  bound  on  BQ.  Thus  we  can  expect  from  our  results  alone:  30.86  s  BQ  s  31.03 
and  3.90  £  B  1  £  4.11.  These  results  have  been  published. 

The  second  pressure  derivative  of  the  bulk  modulus,  Bo",  exhibits  some 
dispersion  and  uncertainty  around  the  value  of  -0.1 (kbar)"1 ,  depending  on  the 
equation  of  state  used,  as  shown  in  Table  I.  This  is  expected  from  the  fact 
that  Bo"  involves  a  third  derivative  taken  from  the  pressure-volume  data.  The 
meaningful  determination  of  BQ"  requires  not  only  an  extremely  accurate 
volume  measurement  with  an  error  aV/Vq  £  10"6,  as  pointed  out  by  Macdonald  and 
Powell  [J.  R.  Macdonald  and  D.  R.  Powell,  J.  Res.  NBS  25A,  441  (1971)],  but 
also  an  equal  accuracy  in  the  pressure  measurements,  as  indicated  by  Kim, 
Chhabildas  and  Ruoff  [K.  Y.  Kim,  L.  C.  Chhabildas  and  A.  L.  Ruoff,  J.  Appl. 

Phys.  47,  2862  (1976)].  The  former  condition  is  certainly  fulfilled  in  our 
volume  measurement,  but  the  latter  is  extremely  difficult  to  meet  in  the 
laboratory. 

The  Grover,  Getting,  Kennedy  equation  [R.  Grover,  I.  C.  Getting  and 
G.  C.  Kennedy,  Phys.  Rev.  B  567  (1973)]  has  been  proposed  as  a  two  parameter 
equation  involving  BQ  and  80'.  Higher  derivatives  of  the  bulk  modulus  are  not 
zero  but  can  be  directly  determined  from  BQ  and  BQ';  we  find 
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B  '* 


which  gives  B0"  s  -0.54.  This  is  too  large  by  a  factor  of  four  to  seven; 
thus  GGK  is  incorrect  on  a  physical  basis  (just  as  ME-j  is);  nonetheless  both 
equations  have  some  practical  use. 

Instead  of  trying  to  meet  the  latter  condition,  the  accuracy  of  BQ"  can 
be  improved,  as  far  as  pressure  is  concerned  by  extending  the  applied  pressure 
range  up  to  pressure  equal  to  Bq  (kbars)  of  a  specimen  or  even  beyond.  A 
potassium  crystal  has  BQ  s  30  kbar.  A  30  kbar  vessel  was  constructed  to 
carry  our  length  measurements  out  to  this  pressure. 

The  overall  view  of  a  30  kbar  vessel  is  displayed  in  Figure  1.  The  two 
stage  pressure  vessel  design  has  been  adopted  to  generate  a  maximum  pressure 
of  30  kbar  in  the  bore  of  the  inner  vessel.  These  vessels  were  made  of  a 
non-magnetic  titanium  alloy,  the  product  of  Reactive  Metal  Corporation, 

RMI  38-6-44.  The  commonly  used  manganin  wire  pressure  sensor  and  a  new  type 
of  pressure  sensor  built  with  a  chip  size  stable  Zener  diode  have  been  construc¬ 
ted  and  were  tested  in  order  to  accurately  calibrate  and  measure  pressure  up 
to  30  kbar.  Improved  feed-through  of  the  electrical  leads  into  30  kbar 
vessel  has  been  devised  to  electrically  detect  the  change  under  pressure. 

Unfortunately,  the  contract  was  not  renewed  and  came  to  an  end  just  as 
the  vessel  was  being  installed  and  tested. 
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THEORETICAL 


I :  Equilibrium  and  Transport  Properties 

The  major  emphasis  of  the  early  part  of  the  theoretical  program  centered 
on  the  equation  of  state  and  transport  properties  of  both  potassium  and 
aluminum.  The  objective  of  the  program  had  been  to  study  the  suitability  of 
the  simple  metals  as  possible  monitors  of  pressure  in  ultra-pressure 
environments.  Under  such  conditions  it  is  likely  that  resistivity  can  be 
measured,  and  for  this  reason  we  concentrated  our  attention  on  electrical 
transport  properties.  These  are  usually  calculated  as  functions  of  volume, 
whereas  the  quantity  of  most  practical  significance  is  the  resistivity  (or 
resistance)  as  a  function  of  pressure.  It  is  therefore  necessary  to  know  the 
equation  of  state  of  the  metal,  particularly  under  conditions  of  high 
compression.  As  described  in  the  original  proposal,  we  pursued  this  aim  on  a 
metal  with  relatively  low  bulk  modulus;  potassium  has  a  zero  pressure  bulk 
modulus  Bq  of  30  Kbar  and  we  examined  the  theoretical  equation  of  state 
at  pressures  far  in  excess  of  30  Kbar.  The  ability  of  the  current  theory  of 
the  metallic  state  to  describe  the  equation  of  state  of  a  metal  for  pressure 
greater  than  BQ  in  a  soft  metal  (such  as  potassium)  was  taken  as  evidence 
that  the  same  theory  would  be  successful  in  describing  a  simple  metal,  such 
as  aluminum,  with  a  much  higher  bulk  modulus. 

These  two  metals  do  differ,  however,  in  an  important  respect.  The  ions 
of  aluminum  are  compact  and  very  tightly  bound.  In  comparison,  the  ions  of 
potassium  are  more  extended  and  less  tightly  bound  electronically.  As  a 
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consequence,  core-polarization  effects  in  potassium  cannot  be  neglected  in  an 
accurate  determination  of  the  equation  of  state.  Thus  in  our  detailed  study 
(1)  we  re-examined  the  theory  of  the  ground  state  energy  of  simple  metals  and 
included  the  electron-electron  interaction  effects  between  localized  core¬ 
states.  This  meant  modifying  the  method  of  structural  expansions  (2,3)  to 
include  the  effects  of  core-polarization  fluctuations.  It  is  well  known  that 
in  insulating  systems  such  fluctuations  lead  to  Van  der  Waals  interactions  at 
long  range.  In  a  metal,  however,  these  interactions  are  screened  by  the 
valence  electrons,  and  since  the  processes  involved  are  dynamic,  the  screen¬ 
ing  requires  a  knowledge  of  the  frequency  dependent  dielectric  function  of 
the  interaction  electron  gas.  That  even  the  screened  Van  der  Waals  interac¬ 
tions  in  potassium  are  appreciable  at  high  pressures  can  be  attributed  to  the 
fact  that  the  pol arizabi 1 ity  of  the  potassium  ion  is  rather  substantial.  The 
results  (1,4)  showed  that  although  it  is  reasonable  to  ignore  core¬ 
polarization  at  low  pressure  (P  <  10  Kbar)  they  become  progressively  more 
important  at  higher  pressures;  for  example,  they  contribute  about  15  Kbar 
(out  of  approximately  60)  at  40%  compression  (5). 

In  addition  to  the  question  of  the  equation  of  state  we  also  examined 
quite  closely  the  behavior  of  the  high  temperature  resistivity  of  aluminum 
under  pressure  (6).  To  minimize  errors  as  much  as  possible,  we  considered 
the  scaled  resistivity,  i.e.,  p(P)/p(P=0),  by  determining  p  from  a  standard 
variational  treatment  of  the  Boltzmann  equation  in  which  the  necessary  elec¬ 
tron  levels  and  distortion  of  the  actual  Fermi  surface  were  obtained  from  a 
replicated  two  plane-wave  model.  The  scaled  resistivity  that  resulted 
displayed  a  very  clear  minimum  (6)  at  a  pressure  of  about  25  GPa.  We  were 
able  to  understand  this  behavior  in  terms  of  the  increasing  importance  of 
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distortions  to  the  Fermi  surface  as  pressure  itself  increased.  In  this 
respect,  the  results  for  the  crystal  are  very  different  from  those  we 
obtained  (7)  for  1 i quid  aluminum  which  possesses  a  spherical  Fermi  surface. 
And  since  the  distortions  to  the  Fermi  surface  in  the  crystal  reflect  the 
changes  in  pseudopotential  with  density,  the  behavior  among  different 
crystalline  simple  metals  is  also  quite  varied  (we  predict  lead  to  have  a 
p(P)  rather  different  from  that  found  in  A1 ) . 

II:  Gruneisen  Parameter  and  Dynamically  Compressed  Metals 

A  major  part  of  the  theory  component  concerned  the  direct  study  of  the 
Gruneisen  parameter,  the  quantity  required  in  reducing  the  shock  data  to  iso¬ 
thermal  conditions.  To  obtain  this  quantity  with  precision  we  required,  as 
functions  of  temperature  and  volume,  the  bulk  modulus,  the  heat  capacity,  and 
the  coefficient  of  volume  expansion. 

As  described  above  we  made  some  advances  in  understanding  the  equation 
of  state  of  potassium  at  high  density,  our  major  interest  was  directed  toward 
the  T  =  0  isotherm.  Subsequently  we  examined  the  contribution  of  the  phonons 
to  the  free  energy  of  a  model  solid  using  the  method  of  self-consistent 
phonon  theory  (7).  This  approach  allows  us  to  build  in  the  effects  of  the 
most  important  anharmonicities,  and  while  these  are  reasonably  small  effects 
at  low  temperatures,  they  become  increasingly  important  at  high  temperatures 
(for  instance,  in  the  conditions  prevailing  in  shock  experiments).  Another 
method  we  developed  (and  are  currently  improving)  gives  a  different  approach 
to  the  problem.  In  this  new  method  (9)  the  thermodynamic  properties  of 
simple  solids  at  high  pressures  are  calculated  from  liquid  state  thermo¬ 
dynamics  by  treating  the  solid  as  a  super-cooled  liquid.  In  particular. 
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we  are  able  to  obtain  the  energy  and  pressure  at  arbitrary  temperatures  and 
densities  via  expressions  that  are  reasonably  simple  functionals  of  the  pair 
potential.  As  ar  application  of  this  principle,  we  have  started  to  analyze 
shock-wave  data  via  liquid  state  thermodynamics.  For  example,  we  can  use 
shock-wave  data  directly  as  a  test  of  a  given  pair  potential  without  re sort¬ 
ing  to  the  usual  techniques  for  "reduction  of  the  Hugoniot."  Alternatively, 
the  method  can  be  used  to  perform  the  reduction  of  a  shock  Hugoniot  by 
selecting  from  a  general  space  of  pair  potentials  the  one  that  best  fits  the 
known  equilibrium  thermodynamic  quantities  and  the  Hugoniot  curve.  The 
resulting  pair  interaction  can  subsequently  be  used  to  generate  zero  tempera¬ 
ture  data  thereby  also  performing  a  "reduction." 

Ill:  Interatomic  Forces  in  Transition  Metals 

We  extended  our  early  work  supported  by  the  grant  on  simple  metals  to 
transition  metals  to  the  extent  that  their  electronic  structure  could  be 
reasonably  well  described  by  the  tight  binding  approximation.  In  particular, 
we  derived  a  new  method  (10)  for  calculating  the  interatomic  potentials  that 
can  be  attributed  to  band  broadening  contributions  to  the  metallic  cohesive 
energy.  Our  method  was  based  on  the  direct  evaluation  in  r-space  of  the 
moments  of  the  electronic  density  of  states,  but  projected  on  a  particular 
atom.  We  applied  it  to  tight-binding  model  Hamiltonians  appropriate  to 
either  one-  or  five-states  per  atom.  In  the  latter  case,  which  is  appro¬ 
priate  to  the  d-orbitals  in  transition  metals,  the  resulting  interatomic 
potentials  have  spatial  forms  rather  characteristic  of  fluctuating  multipole 
forces.  The  band  broadening  contributions  are  attractive:  the  expected 
short-range  repulsive  contributions  are  obtained  by  expanding  the  energy 
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to  second  order  in  atomic  orbital  overlaps.  The  combination  of  attractive 
and  repulsive  contributions  forms  the  basis  of  our  method  for  calculating  the 
thermodynamic  functions  in  the  transition  metals  and  parallels  the  nearly 
free  electron  approach  we  used  in  the  simple  metals.  So  far,  we  have  applied 
the  method  to  the  study  of  vacancy  formation  energies. 

IV:  Electron  Distributions  Around  Charge  Centers  in  Solids 

Finally,  in  the  course  of  studying  pair  and  multi-center  interactions  in 
solids,  we  were  led  to  consider  the  nature  of  electron  distributions  around 
point  charges  in  many-body  systems.  We  managed  to  prove  a  cusp  theorem  (11) 
which  relates  the  zero-separation  value  and  slope  of  two-particle  positional 
correlation  functions  in  quantum  many-body  systems  with  Coulombic  interac¬ 
tions.  The  theorem  is  independent  of  particle  type  and  symmetry  of  the 
wave-function.  Our  proof  uses  only  the  integral  form  of  the  Schrodinger 
equation  and  continuity  and  long-range  exponential  decay  of  the  wave  func¬ 
tion.  We  have  used  the  theorem  to  derive  a  sum-rule  for  the  electron  gas 
structure  factor,  and  an  exact  statement  about  the  screening  of  point 
charges. 
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Volume  versus  pressure  data  of  1-m-long  polycrystalline  potassium  has  been  obtained  as  a 
function  of  hydrostatic  pressure  up  to  7  kbarsat  28.38  °C.  The  length  change  has  been  measured 
to  an  accuracy  of  less  than  500  A  using  a  Fabry-Perot  type  He-Ne  laser  interferometer  in  a 
temperature-controlled  environment  with  temperature  variations  ±  0.001  °C.  The  isothermal 
bulk  modulus  B„  and  its  pressure  derivative  B'0  at  atmospheric  pressure  and  28.58  °C  are 
Ba  =  31.01  +  0.02  kbars  and  B  „  =  3.91  ±  0.01,  when  the  first-order  Mumaghan  equation  is 
used;  and  B0  =  30.84  +  0.02  kbars  and  B'0  =  4. 16  ±  0.02,  when  the  Keane  equation  is  used.  The 
Keane  equation  yields  /?„  =  —  0. 145  kbar' 1  in  the  applied  pressure  range.  Various  two-  and 
three-parameter  equations  of  state  have  been  used  to  fit  the  measured  pressure-volume  data.  The 
Keane  and  Birch  equations  represent  the  data  most  closely  when  these  equations  are  extrapolated 
into  a  higher-pressure  region. 

PACS  numbers-.  64.30.  -t- 1,  64.10.  +  h,  62.50.  +  p 


INTRODUCTION 

Potassium  and  the  other  alkali  metals  have  drawn  con¬ 
siderable  theoretical  interest.  The  measurement  of  their 
thermodynamic  properties  is  helpful  in  understanding  the 
cohesion  of  simple  metals.  However,  their  mechanical  soft¬ 
ness  and  strong  chemical  reactivity  have  made  these  mea¬ 
surements  very  unattractive,  since  extreme  care  must  be  tak¬ 
en  in  handling  these  materials. 

Earlier  measurements  of  volume  compression  of  alkali 
metals  were  carried  out  by  Bridgman 1-4  up  to  100  kbars.  The 
pressure  scale  hp  used  in  a  higher-pressure  region  turned  out 
to  deviate  appreciably  from  the  present  pressure  scale.  Swen¬ 
son  and  his  coworkers5,6  measured  volume  change  up  to  20 
kbars  over  a  temperature  range  down  to  liquid-helium  tem¬ 
perature.  Kennedy  and  his  collaborators7  have  measured  the 
pressure- volume  isotherms  up  to  45  kbars.  Both  of  the  above 
authors  used  the  piston-cylinder  device.  But  the  quasihy¬ 
drostatic  media  and  corrections  due  to  the  piston  compres¬ 
sion  and  friction  contribute  to  errors  in  their  measurements. 
Our  present  measurements  are  under  purely  hydrostatic 
conditions. 

Let  V  denote  the  volume  of  a  specimen  under  pressure  P 
and  at  some  constant  temperature  T.  Then  an  isothermal 
bulk  modulus  B  is  defined  and  expanded  in  terms  of  pressure 
as 


B--vm 

\dP)r 

=  B„  +  b;,{P  -  P„)  +  \B"{P  -  P„)2  +  -  , 


(1) 


where 


and 

are  evaluated  at  atmospheric  pressure  P„. 


Although  most  of  those  measurements  made  by  the 
above  workers  are  quite  adequate  to  describe  thermodynam¬ 
ic  behaviors  at  high  pressures,  they  are  not  precise  enough  to 
derive  reliable  higher-order  pressure  derivatives  of  bulk  mo¬ 
dulus,  such  as  B„  or  higher.  The  potassium  crystal  is  an 
extremely  soft  solid  with  B„  s*  30  kbars  and  B is  expected  to 
contribute  appreciably  to  the  value  of  B  at  higher  pressures. 

Accurate  measurements  of  adiabatic  bulk  modulus  by 
the  ultrasonic  method  have  been  made  by  Smith  and  his 
coworkers."  The  ultrasonic  method  yields  adiabatic  data 
which,  if  isothermal  results  are  desired,  must  be  converted 
by  resorting  to  Overton's  relation.*'  The  Overton  formulas 
involve  a  number  of  parameters,  all  of  which  are  not  accu¬ 
rately  determined  in  the  desired  pressure  range. 

Given  our  precise  length  change  measurements  with 
resolution  better  than  500  A  for  a  1-m-long  potassium  speci¬ 
men.  5„.  B'„,  and  B  <"  have  been  determined  by  statistically 
fitting  the  data  to  several  isothermal  equations  of  state.  A 
detailed  treatment  and  discriminations  between  these  equa¬ 
tions  have  been  described  by  Macdonald  and  his  cowork¬ 
ers. 11,1 1  Macdonald  and  Powell1 1  have  shown  that  it  is  statis¬ 
tically  insignificnat  to  discriminate  the  equations  of  state  in  a 


FIG.  1.  Schematic  of  length  measurement  system.  1:  laser  path  vacuum 
bellows.  2:  laser  beam.  }:  interferometer  mirrors,  4.  coupling  plates.  5: 
LVDT's,  6:  pressure  vessel.  7 :  specimen.  8:  magnetic  cores,  9:  stainless-steel 
tube.  10:  Pyrex  glass  tube. 


24S 


J.  Appl.  Phys.  52(1),  January  1981 


002 1  -8979/8 1/0245-0550 1 . 1 0 


©  1981  American  Institute  o(  Physics 


245 


FIG.  2.  Extrapolation  of  pressure-volume  data  for  potassium.  Bridgman's 
pressure  scale  correction  above  35  kbars,  40—38  kbars,  60—56  kbars, 
80—70  kbars,  and  100—84  kbars.  The  shock  data  are  for  liquid  metal. 


low-pressure  region  (P/B0< 0.01)  unless  the  experimental 
data  are  very  accurate  within  an  error  of  A  V/V0<  10".  The 
precision  in  our  length  change  measurement  is  better  than 
5x10-*  and  P/Bo<0.25. 

EXPERIMENTAL  METHOD 

A  detailed  description  of  the  experimental  setup  has 
been  provided  by  Lincoln  and  Ruoff. 1 2  As  shown  in  Fig.  1 ,  a 
polycrystalline  potassium  rod  which  is  approximately  1  m 


long  and  7.9  mm  in  diameter  was  inserted  in  a  tube,  which  in 
turn  is  held  inside  the  12.7-mm-diam  bore  of  the  nonmagne¬ 
tic  titanium  alloy  pressure  vessel.  The  central  portion  of  the 
tube  is  304  stainless  steel  and  the  remaining  part  of  the  tube, 
where  the  magnetic  cores  of  the  linear  variable  differential 
transformer  (LVDT)  are  slightly  spring-loaded  against  both 
ends  of  the  potassium  specimen,  is  a  Pyrex  glass  tube  which 
is  used  to  avoid  any  spurious  magnetic  effects  on  the  H- 
shaped  magnetic  cores. 

Potassium  was  obtained  from  Mine  Safety  Appliance 
Corporation  and  is  99.95%  pure.  This  was  extruded  into  a 
7.9-mm  die  to  the  desired  length.  It  is  then  zone  refined  sev¬ 
eral  times  to  enhance  purity  and  to  anneal  out  defects  and 
some  dislocations.  A  more  elaborate  ball  bushing  arrange¬ 
ment  for  holding  the  specimen  to  minimize  friction  was  pre¬ 
viously  used.12-14  This  was  not  adopted  in  this  experiment, 
because  the  small  difference  in  density  between  potassium 
and  the  pressurizing  medium  (liquid  hexane)  keeps  friction 
sufficiently  small.  Dehydrated  liquid  hexane  is  used  as  the 
pressure  fluid.  Pressure  was  measured  by  using  a  manganin 
wire  pressure  sensor.  The  calibrating  technique  of  the  man¬ 
ganin  gauge  is  described  in  detail  elsewhere.13 

The  length  measurements  were  made  at  an  interval  of 
300  bars  up  to  a  maximum  pressure  of  7  kbars.  Inasmuch  as 
the  linear  compressibility  of  polycrystalline  potassium  is  iso¬ 
tropic,  volume  measurements  are  obtained,  using  the  rela¬ 
tion  V/V0  =  (///0)3.  Resolution  in  the  length  measurements 
of  the  He-Ne  laser  interferometry  is  better  than  300  A  in  the 
temperature  stabilized  environment  with  thermal  noise  of 
0.001  *C. 

The  LVDT  electronics,  laser  interferometer,  coupling 
plates  and  their  translation  motion,  and  temperature  con¬ 
trol,  etc.,  have  been  described  in  detail  in  the  literature.12 

DATA  ANALYSIS 

For  notational  convenience,  we  shall  introduce  the  fol¬ 
lowing  abbreviations: p  =  P~P0,Tj  —  B„,  V  =  BqB £, 
z  =  P/Bot  and  x  =  Vg/V.  Note  that  17,  IP,  r,  and  x  are  all 
dimensionless  quantities.  Po  is  assumed  to  be  1  bar  in  this 


TABLE  I.  Various  equations  of  state. 


Equation 

Acronym 

Form 

First-order 

ME, 

*  =  l| 

Mumaghan 

*  •  2(jc<t?,'2'^),/I  -  l)AtV2  - 

Second -order 

ME, 

Mumaghan 

( 

Keane 

KE 

*  =  V1  +  *)J] 

(  -  »7J  <  46  <  0) 

Xi*1*’*"',-  l|  —  [*/(,» +  *)]!„* 

First-order 

BE, 

*  =  +  ((  V  -  4X*J/’  -  I)] 

Birch 

Second -order 

BE, 

1  =  tf*”’  -  a5'5)!  1  +  j(  r,  -  «X  x21'  -  1) 

Birch 

+  AH«  +  *><V-7)-l-9*ltJrJ,J-l|I| 

Grover 

OGKE 

Br  =  B„r espial  -x  ’)l 

Getting, 

Kennedy 
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TABLE  II.  Bulk  modulus  and  its  pressure  derivatives  for  potassium  at  28.S8  *C. 


Equation 

Bo 

Bo 

Used 

(kbar) 

B'o 

(kbar*1) 

ME, 

31.01  ±0.02 

3.905  ±  0.009 

0 

ME. 

30.87  ±  0.01 

4.094  ±  0.018 

-  0.0747  ±  0.0068 

GGKE 

30.81  ±0.01 

4.184  ±0.004 

-  0.568*  ±  0.001 

BE, 

30.85  ±  0.01 

4.144  ±0.003 

-0.131  *±0.002 

BE, 

30.85  ±  0.02 

4.151  ±0.021 

—  0.135  ±0.010 

KE 

30.84  ±  0.02 

4.159  ±0.024 

-0.145  ±0.014 

‘Obtained  from  0O  and  B'0. 


paper.  Several  of  the  two-  and  three-parameter  phenomolo- 
gical  equations  described  in  detail  elsewhere10, 16  are  used  to 
fit  the  pressure-volume  data  and  are  listed  in  Table  I. 

Table  II  shows  the  values  of  the  bulk  modulus  and  its 
pressure  derivatives  obtained  by  fitting  the  measured  pres¬ 
sure  volume  data  into  the  equations  of  state  listed  in  Table  I. 
ME,  is  the  only  equation  of  state  in  the  list  for  which  B„  is 
zero.  The  expressions  for  BE,  and  GGKE  do  not  explicitly 
contain  the  parameter  B  „ ;  however,  this  does  not  imply  that 
B  o  is  identically  zero.  The  appropriate  expression  for  B  „  in 
the  case  of  BE,  can  be  obtained  by  setting  the  third  term  in 
the  second-order  Birch  equation  BE,  equal  to  zero.  The  re¬ 
sult  is 


9B„B;  =  -  [143  +  9fi'(£'  -7)].  (2) 

Successively  differentiating  the  GGKE  equation  with  re¬ 
spect  to  pressure  and  taking  the  limit  as  /* — ►/*,„  one  obtains 

W  =  -fl,'2.  (3) 


GGKE,  along  with  several  other  equations  of  state,  fits  the 
pressure-volume  data  well.  However,  a  low  value  of  B'a  (30.8 
kbars)  gives  an  unreasonably  high  value  of  B„,  as  expected 
from  Eq.  (3). 

Another  approach  for  the  determination  of  the  bulk 
modulus  at  each  pressure  point  is  to  combine  three  neighbor¬ 
ing  data  points  and  fit  them  into  ME,.  That  is, 

^ _ .  ,4, 

V  Jr.  B,  +  B  '(P  -  P, ) 

^=f  • - ^ _ .  (5, 

.  V  Jr.  B,  +*;(/>-/>,) 


Out  of  the  two  relations  (4)  and  (5),  B,  and  B '  have  been 
numerically  calculated  at  each  pressure  point  (except  P  =  P0 
and  P  =  Pm„ ).  These  bulk  modulus  versus  pressure  data  are 
statistically  fitted,  using  the  formula 


The  values  of  B0  and  B  for  N  =  1 , 2,  and  3  are  displayed  in 
Table  III. 


As  expected,  ME,  in  Table  II  and  N  =  1  in  Table  III, 
yield  essentially  the  same  result  for  B0  and  B'0,  and  ME2  in 
Table  II  and  N  =  2  in  Table  III  also  show  essentially  the 
same  result  for  B0,B'0,  and  B  within  experimental  errors, 
respectively.  This  provides  a  self-consistency  check  for  ME, 
and  ME,.  It  is  interesting  to  notice  that  N  =  3  inTable  III 
provides  almost  identical  results  for  B0,B  „ ,  and  B  „  to  those 
which  BE,  and  KE  yield  (Table  II).  Note  that  the  negrect  of 
B„  in  ME,  leads  toaB„  value  which  is  too  small. 

Table  IV  shows  the  data  of  bulk  modulus  and  its  pres¬ 
sure  derivatives  measured  by  other  authors  together  with  the 
present  work  for  the  sake  of  comparison.  The  agreement  in 
B  o  is  only  adequate  when  compared  with  Bridgman  and 
Kennedy's  static  compression  data,  but  B0  and  B  g  are  in 
good  agreement  with  Swenson’s  static  compression  data  and 
Smith's  values  obtained  by  the  ultrasonic  technique. 

DISCUSSION 

As  can  be  seen  in  Tables  II  and  IV,  the  two-  and  three- 
parameters  equations  yield  virtually  identical  results  for  B„, 
and  values  of  B  o  obtained  from  GGKE,  BE,,  BE2,  and  KE 
are  also  in  good  agreement.  Although  any  sign  of  curvature 
is  hardly  noticeable  in  the  ultrasonic  mesurementsof B  (P )  vs 
P  for  potassium,"  the  three-parameter  equations  (ME2,  BE,, 
KE)  show  unmistakably  a  negative  value  of  Bp  a  —  0.1 
kbar  '.  Anderson17  used  Swenson’s  value5,6  for  B0  and  B'0 


TABLE  III.  Values  of  0„  and  B‘,"‘  for  ,V  =  1,2,  and  3 


B„ 

N  (kbar) 

0' 

B", 

(kbar  ') 

Bo1 

(kbar'J| 

1  31.06  ±0,05 

3.88  ±  002 

0 

0 

2  30.86  ±  0.04 

4.09  ±  0.03 

-  0.073  ±  0.010 

0 

3*  30.81  ±0.05 

4.18  t  0.08 

-0.145  ±0063 

0.024  ±0.021 

‘Because  of  the  large  standard  deviation  associated  with  0.7  '  and  the  much  larger  standard  deviation  associated  with  0,7  than  with  the  V  •  2  case,  the 
numbers  in  this  row  should  not  be  given  great  significance. 
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TABLE  IV.  Isothermal  bulk  modulus  and  its  pressure  derivatives  of  potassium  (equation  used:  ME,,  MEJt  ME,). 


Temperature 

Bo 

B~n 

B o' 

Source 

ra 

(kbar) 

Be 

(kbar1) 

(kbar  j 

Bridgman4 

Swenson’ 

Room 

33.25 

3.43 

25* 

30.8 

3.85 

Smith4  * 

Room 

30.9 

3.98 

Kennedy’ 

Room 

34.0 

2.99 

Room 

31.2 

3.65 

-  0.039 

This  work b 

28.58 

31.01 

3.91 

28.58 

30.87 

4.09 

-0.075 

28.58 

30.81 

4.18 

-0.145 

0.024 

‘Corrected  to  an  isothermal  value  from  ultrasonic  data. 
'’See  comment  in  Table  III  regarding  last  row  of  data. 


taken  at  low  pressures  on  the  solid  and  chose  B  7  which  can 
best  tit  the  shock  wave  data'8  (which  is  for  the  liquid)  taken  at 
higher  pressures,  using  the  Keane  equation.  He  obtained  the 
values  of  B  7  =  —  0.05 1  kbar  “ 1  for  potassium  and  B  7 
=  —  0.031  kbar- 1  for  sodium.  Similar  procedures  have 
been  used  by  Fritz  and  Thurson by  using  a  different  type  of 
extrapolation  formula.  Their  values  of  Bq  are  —  0.018 
kbar  ~ 1  for  potassium  and  0.024  kbar  ~ 1  for  sodium.  Consid¬ 
ering  that  this  involves  mixing  shock  data  on  the  liquid  with 
low-pressure  isothermal  data  on  the  solids,  these  values  of 
5,7  should  not  be  taken  seriously. 

The  value  of  5,7  shows  very  sensitive  dependence  not 
only  on  the  accuracy  of  the  V /V0  measurement,  as  pointed 
out  by  MacDonald  and  Powell,1 1  but  also  on  the  functional 
form  of  the  variation  of  pressure  with  the  change  in  resis¬ 
tance  per  unit  resistance  AR  //?„  of  the  manganin  gauge,  as 
indicated  in  the  data  of  LiF  and  NaCl.1  ,M  Our  pressure 
measurements  are  known  to  a  precision  of  1 X  10  ~4.  Unless 
the  pressure  is  measured  to  an  extremely  high  precision  of 
1 X  10  A,  the  error  in  57  is  going  to  remain  relatively  large 
unless  the  measurements  are  made  to  substantially  higher 
pressure. 

The  extrapolation  of  the  pressure-volume  data  into  a 
higher-pressure  range  is  displayed  in  Fig.  2  in  order  to  see  the 
discrepancies  between  the  equations  of  state,  despite  the 
questionable  validity  of  extrapolation.  Bridgman's  pressure 
scale  above  35  kbars  was  corrected  by  using  the  recently 
reported  values  of  pressure  transition  points  of  various  ele¬ 
ments.'"  His  values  of  thallium  II-III  at  39  kbars,  barium  II- 
III  at  59  kbars,  and  bismuth  III-V  at  89  kbars,  moved  down 
to  37, 21  55,22  and  74-78  kbars,2'24  respectively.  GGKE  has 
not  been  shown  because  of  its  unrealistically  high  value  of 
5,7.  As  shown  in  Fig.  2,  ME,  and  ME,  begin  to  deviate 
appreciably  from  other  equations  of  state  at  30  kbars.  The 
abnormal  behavior  of  ME,  in  the  higher-pressure  region  is 
due  to  the  dominant  role  of  a  5,7  term  in  the  determination 
of5{5).  KE,  BE,,  and  BE,  are  hardly  distinguishable  in  the 
entire  pressure  range  (5/5, , a  10)  considered. 

For  the  improvement  of  equation-of-state  data  it  would 
be  desirable  to  simultaneously  measure  pressure-volume 
data  by  the  present  apparatus  and  5,  (5 )  by  an  ultrasonic 
technique  for  the  same  specimen.  Accurate  measurements 


over  an  extended  pressure  range  will  certainly  contribute  to 
improving  data  in  the  determination  of  57  and  57 

CONCLUSIONS 

The  value  of  5„  for  potassium  at  28.58  *C  is 
30.85  ±  0.05  kbars. 

The  value  of57  for  potassium  at  28.58  °C  is4.1  +0.1. 

5  7  >s  found  to  be  5  7  at  —0.1  kbar  - 1  and  hence  makes 
a  major  contribution  to  the  value  of  5  at  high  pressures.  It 
appears  that  the  Taylor  series  expansion  of  5  (5 )  about  atmo¬ 
spheric  pressure  is  a  slowly  converging  series.  This  is  consis¬ 
tent  with  the  general  conclusions  of  Davison  and  Graham25 
that  fourth-order  elastic  constants  contribute  significantly 
to  the  stress-strain  response  even  at  strains  of  only  a  few  per 
cent. 
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The  scaled  room-temperature  resistivity  [p(p)/p(p  =  0))  of  crystalline  aluminum  is  calculated  as  a  function 
of  pressure  p.  Initially  the  resistivity  is  determined  as  a  function  of  volume  from  the  standard  variational 
treatment  in  which  the  required  electron  levels  and  distortions  to  the  Fermi  surface  are  described  in  a  two- 
ptane-wave  model.  To  obtain  the  resistivity  as  a  function  of  pressure,  the  results  of  this  calculation  are 
combined  with  a  previously  computed  equation  of  state  aluminum.  The  calculated  scaled  resistivity  then 
shows  a  minimum  at  a  pressure  of  about  25  GPa.  This  minimum  is  largely  attributable  to  the  increasing 
importance  of  distortions  of  the  actual  Fermi  surface  as  pressure  increases. 


I.  INTRODUCTION 

The  electronic  structure  of  aluminum  is  rela¬ 
tively  simple.  Its  bands  have  a  largely  free-elec¬ 
tron-like  character  and  can  be  interpolated  quite 
accurately  by  a  spatially  local  pseudopotential. 

On  the  other  hand  its  Fermi  surface  is  a  complex 
multiply  connected  object  that  is  sensitive  to  the 
choice  of  the  pseudopoteiuial  components  used  to 
interpolate  the  band  structure.  Since  the  trans¬ 
port  coefficients,  and  in  particular  the  resistivity, 
are  related  to  integrals  over  the  Fermi  surface, 
one  might  expect  this  sensitivity  to  become  ap¬ 
parent  if,  as  through  the  application  of  pressure, 
the  pseudopotential  coefficients  are  altered.  As 
we  shall  see  below,  this  is  indeed  partly  the  case, 
though  in  crystalline  aluminum  the  pressure  de¬ 
pendence  of  resistivity  turns  out  to  be  an  aggre¬ 
gate  of  some  partially  compensating  effects.  This 
compensation  is  very  much  a  property  of  the  me¬ 
tal  itself  (in  Pb,  for  example,  the  effects  we  dis¬ 
cuss  should  be  more  prominent)  and  also  of  its 
state  (in  solid  aluminum  the  effects  are  far  more 
noticeable  than  in  the  liquid  state1). 

We  are  concerned  in  this  paper  (as  in  Ref.  1) 
with  the  scaled  room-temperature  resistivity 
p(p)  'pip  =0)  of  crystalline  Al,  at  a  pressure p. 

The  natural  quantity  to  calculate  is  the  ratio 
p(0  pig,  where  Vis  the  volume  of  a  sample 
at  pressure  p  (and  1’0  its  value  at  p  =0).  The  start¬ 
ing  point  of  this  calculation  is  the  well-known  var¬ 
iational  expression*  for  a  bound  on  p(V),  as  de¬ 
scribed  in  Sec.  n.  For  metals  with  complicated 
Fermi  surfaces,  the  necessary  computations  gen¬ 
erally  require  numerical  procedures  of  matching 
complexity,  even  for  relatively  simple  choices  of 
the  variational  trial  function.  The  result  for  p(V) 
can  certainly  be  expected  to  depend  on  this  choice 
as  well  as  on  approximations  made  necessary  for 
wholly  numerical  reasons.  Much  of  the  consequent 
uncertainty  can,  however,  be  reduced  by  focusing 
attention  on  the  scaled  quantity  p(0/p(t'0)  and,  as 


discussed  in  Sec.  II,  it  is  largely  for  this  reason 
that  we  find  it  convenient  to  use  the  simplest  form 
of  trial  function.  Within  this  approximation  it 
still  remains  to  determine  the  behavior  of  the 
pseudopotential  vt,  the  band  structure,  and  the 
Fermi  surface  as  functions  of  volume.  This  is  de¬ 
scribed  in  Sec.  HI.  The  piecewise  two-plane-wave 
approximation  that  we  use  in  evaluating  the  varia¬ 
tional  integrals  is  described  in  Sec.  IV.  The  ma¬ 
trix  elements  appearing  in  the  integrand  also  re¬ 
quire  the  phonon  frequencies  and  their  volume  de¬ 
pendences,  a  question  that  is  taken  up  both  in 
Secs.  IV  and  V.  The  additional  approximations 
we  make  in  order  to  complete  the  numerical  pro¬ 
cedures  are  described  in  more  detail  in  Sec.  IV. 
They  involve  certain  simplifications  in  the  Fermi- 
surface  geometry  and  in  the  description  of  the 
electronic  levels  associated  with  that  geometry. 
The  results  are  discussed  in  Sec.  VI. 

The  calculations  we  report  could  be  performed 
in  principle  for  all  simple  metals.  We  have  se¬ 
lected  aluminum  for  the  reasons  given  earlier,’ 
namely  that  its  high  electron  density  and  small  ion 
core  imply  an  ability  to  sustain  a  high  pressure 
without  core  contact.  In  addition,  the  equation  of 
state  of  Al  has  been  calculated  to  pressures  in 
excess  of  300  GPa  (3  Mbar).  This  information  al¬ 
lows  us  to  convert  from  pOVpO'o) to  p(/>)/p(/>  =0) 
and  hence  arrive  at  the  curves  described  in 
Sec.  V. 

II.  RESISTIVITY  OF  SIMPLE  METALS 

For  a  metal  of  valence  Z,  the  standard  varia¬ 
tional  reduction  of  the  Boltzmann  equation  yields 
for  the  resistivity  at  high  temperatures  the  ex¬ 
pression2 


U’.v'r  !;V't,)4>(.v)r  2PSf.V 1  ’ 
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l€x(%,  -  x,)-(x,  |Vit>  IX,)  i* 
sinh^l^w^Xj-x,)] 


Here  4>  is  (he  trial  function  used  in  obtaining  the 
bound  for  p,  kr  is  the  magnitude  of  the  Fermi 
wave  vector,  vr  is  the  Fermi  velocity,  Sr  is  the 
Fermi  energy,  0=  1  /kaT,  M  is  the  mass  of  an  ion, 
and  wx  and  Sx  refer  to  the  frequency  and  polariza¬ 
tion  of  phonons  with  reduced  wave  vector  x.  We 
use  scaled  wave  vectors  x  =  C/2 and  a  scaled 
pseudopotential  w  =  where  is  the 

known  long-wavelength  limit  of  the  pseudopoten¬ 
tial  form  factor.  The  quantity  a^fi/e2  is  the  atomic 
unit  of  resistivity  and  has  the  value  21.7  pQcm. 

In  the  one-plane-wave  approximation  the  quantity 
/  in  (1)  can  be  written 


x  -  x,)|JS(x,,x) ,  (2) 


where  S(xitx.2)  is  the  one-phonon  structure  which 
takes  the  form 


S(x)  =  saca-xl)=2/35,J  £  (€x(*)-x]J 

x  sinhJ[50/Tu>x(x)]. 

(3) 

The  sum  in  (3)  is  over  the  possible  polarizations 
A,  and  the  integral  is  over  the  actual  Fermi  sur¬ 
face  which  in  the  case  of  the  alkali  metals,  for  ex¬ 
ample,  can  be  very  well  approximated  by  a  sphere. 
For  the  polyvalent  metals,  however,  the  shape  of 
the  Fermi  surface  must  be  taken  into  account 
since  a  substantial  fraction  of  the  free-electron 
Fermi  surface  can  actually  be  lost. 

It  has  been  shown*’5  that  for  the  calculation  of 
high-temperature  resistivity  the  effects  of  aniso¬ 
tropy  in  the  choice  of  trial  functions  are  of  dimin¬ 
ishing  importance.  As  a  consequence  of  this  ob¬ 
servation  we  have  used  a  trial  function  of  the  form 


* (x)<x  v(x)’i  ,  (4) 

(where  z  is  parallel  to  the  electric  field)  rec¬ 
ognizing  that  although  the  calculated  resistivity 
must  necessarily  be  in  excess  of  the  actual  value, 
the  effect  of  the  approximation  on  p(p)/p(0),  as 
noted  earlier,  will  be  much  reduced.  Although  this 
choice  for  4>  greatly  simplifies  the  numerical 
work,  a  generalization  to  more  complex  forms  for 
>f>  is  quite  straightforward. 

Given  (4)  as  the  trial  function  then  for  a  cubic 
system,  (1)  becomes 


where 

m  _  fdsx  y(E)  .g. 

m0pt  sqC/- 

defines  the  optical  effective  mass  which  can  be 
calculated  directly.*  In  (6)  s0  is  the  area  of  the 
free-electron  Fermi  surface  (s0-4vk2F).  Atten¬ 
tion  then  reverts  to  the  remaining  double  surface 
integral  in  /,  which  we  shall  treat  in  a  two-plane- 
wave  approximation.  Were  the  Fermi  surface 
spherical,  a  one-plane-wave  treatment  would  suf¬ 
fice  and  the  identity 

/  dsf  J  dsn. 

*/  *F 

=  J  die  J  dV  6(k-  kF)t>(k'  -  k, )  (7) 

reduces  the  problem  to  a  single  volume  integral 
over  q  =  ic  -  fc'.  With  these  simplifying  features  the 
calculation  of  the  resistivity  in  the  alkali  metals 
is  itself  relatively  simple.  In  the  polyvalent  metals, 
however,  transformation  (7)  will  not  correctly 
treat  transitions  involving  parts  of  the  Fermi  sur¬ 
face  that  depart  from  a  simple  spherical  charac¬ 
ter.  On  the  other  hand  (as  we  shall  see)  there  are 
portions  of  the  actual  Fermi  surface  that  remain 
very  spherical  and  contributions  from  these  can 
be  transformed  into  a  corresponding  volume  inte¬ 
gral.  The  remaining  (nonspherical)  portions 
must,  however,  be  treated  by  evaluating  (1)  direct¬ 
ly. 


III.  PSEUDOPOTENTIAL  AND  FERMI-SURFACE 
GEOMETRY 


In  the  two-plane-wave  approximation  to  be  dis¬ 
cussed  in  Sec.  IV,  the  electronic  energy  $  for  level 
k  is  given  near  zone  planes  by  the  solutions  of 


Sn-S  V,(K) 

>,(K)  Sn.'K-S 


(8) 


where  K  is  the  reciprocal-lattice  vector  under 
consideration  and  (/I*/2m)k*.  V/K)  is  related 
to  the  screened  pseudopotential  ^(K)  by  V,(K) 

=  I7jIi>>(K),  where  fl0  is  the  volume  of  the  primitive 
cell:  by  choice  T’,{0)  =0. 

In  order  to  calculate  the  resistivity  we  need  to 
know  the  form  of  the  pseudopotential  v,{k)  through¬ 
out  the  range  0  <  k  «  2kr  For  aluminum,  a  one- 
parameter  empty-core  potential  serves  adequate¬ 
ly.7  In  £  space  it  can  be  written  (with  x*k/2kr 
and  s  =  2 kfre) 


i>(x)  = 


-r.Ze2 
x1^  (x) 


cos(s.v) . 


(9) 


Here  rt  is  the  parameter  usually  referred  to  as 
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the  “empty-core  radius.”  The  value  of  should 
be  close  to  the  size  of  the  ion  core  on  physical 
grounds  and  is  around  1.12a0  for  aluminum.  For 
the  dielectric  function  e(x),  we  take* 

c  (*)  =  1  +  tfU)[l  -  C(*)] , 
with 

H(x)  =  (»a0ft^s)-1(|+  (10) 

and 

G(x)=x'/(2x2  +  i). 

The  exchange-correlation  parameter  4  is  given  by* 

(  =  l/(l  +  ar,) , 

where  r,  is  the  electron  spacing  radius  defined  by 
n=itZrl,  and  or  as0.025. 

The  screened  pseudopotential  form  factor  for 
aluminum  is  plotted  in  Fig.  1.  In  the  same  figure 
we  also  show  the  corresponding  quantity  for  the 
compressed  metal  with  r'Uo  =  0.7.  On  plotting  w 
as  a  function  of  x  we  see  that  the  node  (where  w 
=  0)  moves  to  lower  reduced  wave  vector  as  the 
metal  is  compressed.  Correspondingly  the  Fourier 
coefficients  V^lll)  and  1^(200)  are  seen  to  in¬ 
crease.  While  pseudopotential  theory  suggests 
that  r,  should  be  energy  dependent,®  the  energy 


FIG.  1.  Pseudopotential  form  factor  of  aluminum  for 
VA\  *  1.0  and  t'/V,  »  0.7.  Screening  by  conduction  elec¬ 
trons  (Eq.  (10)1  the  appropriate  densities  is  taken  into 
account,  the  empty-core  rndlus  re  is  l.llnu.  The 
arrows  mark  the  location  of  the  reciprocal-lattice  vec¬ 
tors  K  •  1,  1)  and  (2,  0,  0). 


dependence  is  expected  to  be  small  since  re  re¬ 
flects  the  size  of  the  ion  core  which  should  not 
change  greatly  with  pressure-induced  changes  in 
its  environment.  Holding  re  fixed  is  tantamount 
to  ignoring  such  small  additional  energy  depen¬ 
dence. 

With  the  "new"  values  of  V,(lll)  and  F,(200) 

(for  the  compressed  metal)  we  can  use  (8)  again 
to  map  out  the  Fermi  surface.  Because  of  the  dif¬ 
ferent  values  of  the  new  Fourier  coefficients,  the 
geometry  of  the  Fermi  surface  of  the  compressed 
metal  can  be  different  from  that  of  the  metal  under 
normal  conditions.  We  take  this  into  considera¬ 
tion. 

IV.  TWO-PLANE-WAVE  APPROXIMATION 

The  evaluation  of  the  double  surface  integral  in 
(1)  has  been  a  major  numerical  obstacle  in  most 
of  the  calculations  of  the  resistivity  of  polyvalent 
metals.  Because  of  the  complex  geometry  of  the 
Fermi  suface  of  these  metals  a  large  number  of 
surface-area  elements  is  necessary  to  charac¬ 
terize  the  surface  accurately.10'12  For  the  high- 
temperature  transport  coefficients  the  problem 
of  the  Fermi-surface  anisotropy  is  less,4'5  and  in 
view  of  this  we  have  chosen  to  carry  out  the  cal¬ 
culations  by  means  of  a  two-plane-wave  approxi¬ 
mation  which  has  also  been  used  by  other  authors 
in  similar  model  calculations.11,14  Although  at 
some  points  on  the  Fermi  surface  three  or  even  . 
four  orthogonalized  plane  waves  are  needed  to  give 
an  adequate  description  of  the  finer  distortions,15 
the  amount  of  surface  requiring  this  more  detailed 
description  is  small  compared  with  the  total  Fer¬ 
mi-surface  area.  Essentially,  our  approximation 
treats  each  of  the  many  Bragg  planes  in  turn  and  . 
calculates  the  contribution  to  the  resistivity  with 
the  electronic  levels  described  by  the  linear  com¬ 
bination  of  plane  waves 

♦g(r)  =  sin$£e,^'r  +cosfige,<^*^,'~.  (11) 

The  total  resistivity  is  then  the  sum  of  the  um- 
klapp  contributions  of  the  individual  reciprocal- 
lattice  vector  (iv)  and  the  normal  contribution. 

Our  numerical  calculations  show  that  the  normal 
contribution  is  only  a  small  part  of  the  total  re¬ 
sistivity  and  for  this  it  is  quite  adequate  to  use 
the  simple  one- plane- wave  treatment. 

Consider  the  umklapp  processes  made  possible 
by  transitions  involving  a  particular  reciprocal- 
lattice  vector  iv.  The  umklapp  processes  are 
those  for  which  the  initial  levels  It,  originate  on 
the  Fermi  surface  and  the  final  levels  It,  end  on  the 
remapped  surfaces,  i.e., 

ft,  -  It,  =  q  -  !v  (12) 

with  q  in  the  first  Brillouin  zone.  Some  umklapp 
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FIG.  2.  Two  possible  umklapp  processes  are  shown 
here  in  thejslngle-Bragg-plane  (BP)  approximation  (kt 
—  k2;  kj  —  kj)-  Note  that  the  Fermi  surface  Is  divided 
up  into  three  parts  for  computational  purposes  as  de¬ 
scribed  in  the  text. 

.  processes  are  pictured  in  Fig.  2  where  we  also 
see  that  the  Fermi  surface  is  split  into  two  sec¬ 
tions  because  of  the  finite  value  of  vj. 

The  parameter  9  describing  the  electron  levels 
in  (11)  is  given  by 

tanSj  =  tj±  sgn(p*)(t/2  + 1)1/2 , 
with 

'?  =  [4S,/|V;(K)|] 

x[Z-(K/2kF)-L(K/2kF)2],  x=k/2kr.  (13) 

The  (±)  in  (13)  signify  the  band  index,  with  (+)  for 
the  second  band  and  (-)  for  the  first. 

The  matrix  element  M  required  in  (1)  is  in  gen¬ 
eral  given  by 

-  ii  x  (Xj  -  x,)  •  <x2  |  Vic  |x,) 

when  the  wave  functions  are  normalized.  Before 
we  consider  this  further  we  will  first  describe  the 
approximations  we  make  concerning  the  phonon 
frequencies  and  polarizations.  With  the  excep¬ 
tion  of  those  along  the  major  symmetry  axes,  the 
phonons  have  polarizations  that  are  neither  pure 
longitudinal  nor  pure  transverse.  However,  the 
large  number  of  symmetry  axes  present  in  the 


Brillouin  zone  of  the  cubic  structures  guarantees 
that  throughout  the  zone  the  phonons  have  polari¬ 
zations  that  are  usually  quite  close  to  being  purely 
longitudinal  or  purely  transverse.  It  thus  seems 
reasonable  to  replace  the  phonon  dispersion  curves 
with  one  longitudinal  branch  and  two  identical 
transverse  branches  with  dependence  on  the  mag¬ 
nitude  of  q  alone,  these  branches  being  appro¬ 
priately  weighted  averages  of  the  frequencies 
along  the  major  symmetry  axes.  Implicit  in  this 
is  the  subsequent  replacement  of  the  zone  by  the 
Debye  sphere.  The  individual  phonon  frequencies 
are  calculated  from  an  eight-shell  axially  sym¬ 
metric  model  with  experimental  force -constants 
obtained  from  experiments  at  300  K.1*  This  pro¬ 
cedure  inevitably  introduces  errors  into  the  final 
result,  but  the  uncertainties  are  not  expected  to 
be  large,  especially  if  the  main  interest  centers 
on  the  way  resistivity  changes  as  a  consequence 
of  the  variation  of  the  frequencies  themselves. 

With  this  approximation,  the  scaled  longitudinal 
and  transverse  matrix  elements  (squared)  are 
then  given  by 

=  {?  ^(q)(S,S2  +  C,C2)  +  7  F,(q  -  K)S2C, 

+  q  F,(q  +  K)S,  C2  +  q  •  K[  l',(q  +  K)S,  C2 

-tyq-KjS^j2!2  (14) 
and 

(2*j-)J(§Sf)2  A/J.fltj.lc,) 

=K2(i-(r^)JJ(v,,(q  +  K)slci 

where  Sf  =  sin0,-  and  Ct  =  cos9,  with  » =  1,2.  Given 
the  trial  function  (4)  the  quantity  [$(F2)  -  *(^i)F  ap¬ 
pearing  in  (1)  is  determined  by 

=  72  +  2q-K(Cj-C2)  +  tf2«^-C2)2.  (15) 

If  the  band  gaps  are  small  (as  in  aluminum)  then 
only  a  small  region  of  the  Fermi  sphere  needs 
to  be  described  in  the  two-plane-wave  formalism, 
the  rest  being  essentially  free-electron-like.  Con¬ 
siderable  computational  effort  can  be  saved  if  we 
divide  up  the  first  band  into  two  parts  (see  Fig.  2), 
/,  and  /t,  in  which  /,  is  the  part  of  the  first  band 
where  a  one-plane-wave  description  would  be  ade¬ 
quate  while  If,  is  the  part  that  requires  the  two- 
plane-wave  description  (as  determined  by  a  suit¬ 
ably  chosen  criterion12).  The  processes  involving 
initial  states  on  /,  and  final  states  on  the  remapped 
part  of  /4  are  essentially  free-electron-like  and 
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will  be  dealt  with  separately.  Contributions  from 
processes  involving  II  and  I„  will  then  be  computed 
by  directly  evaluating  the  double  integral. 

The  double  Fermi-surface  integrals  required  for 
the  evaluation  of  (1)  thus  have  contributions  from 
the  normal  processes  (which  are  small  in  the 
case  of  aluminum),  umklapp  processes  in  the  two- 
plane-wave  single-Bragg-plane  approximation, 
and  umklapp  processes  that  can  be  considered  as 
involving  free-electron-like  levels. 

V.  APPLICATION  TO  COMPRESSED  ALUMINUM 

As  mentioned  above,  the  phonon  frequencies  at 
normal  conditions  (300  ‘X  and  zero  pressure)  are 
generated  by  an  axially  symmetric  force-constant 
model.  The  effects  of  the  volume  change  on  the 
phonon  frequencies  themselves)  can  be  calculated 
in  a  straightforward  manner  from  the  dynamical 
matrices. 16,19  In  the  polyvalent  metals  it  is  well 
known  that  such  a  calculation  is  computationally 
time  consuming  as  a  large  number  of  terms  is 
needed  for  the  dynamical  matrices  to  converge. 

We  choose  instead  to  scale  all  the  frequencies  with 
the  bulk  experimental  Griineisen  parameter  y 
(2.35  for  aluminum20),  i.e., 


>(q,  r)  =  t*J0(<7)^l 


In  fl 6)  u>(<7,  V)  is  the  phonon  frequency  of  the  wave 
vec-  .r  q  at  the  compressed  volume  V  and  w0  is 
the  observed  frequency  at  the  zero-pressure  vol¬ 
ume  K0.  Though  a  crude  approximation  for  the 
changes  in  the  frequencies  themselves,  it  should 
be  satisfactory  for  the  quotient  p(V)/p(V0)  involv¬ 
ing  ratios  of  integrals  accompanying  such  changes. 
A  more  realistic  approach  is  to  take  into  account 
the  changes  in  the  elastic  constants  in  the  evalua¬ 
tion  of  the  changes  in  the  phonon  frequencies.31 

The  numerical  evaluation  of  Eq.  (1)  also  involves 
the  computation  of  the  factor  (m0p,/”^),.  From  Eq. 
(6)  we  see  that  this  factor  is  a  measure  of  the  dis¬ 
tortions  in  the  Fermi  surface.  On  the  other  hand, 
we  note  that  in  evaluating  the  umklapp  contributions 
to  the  double  surface  integral,  the  Fermi  surface 
used  is  not  the  actual  fully  distorted  Fermi  sur¬ 
face.  For  each  particular  Bragg  plane  (at,  say, 

£K)  only  the  distortions  associated  with  a  given 
l'p(K)  are  taken  into  account.  This  amounts  to 
using  a  Fermi  surface  with  an  area  larger  than  the 
actual  one  (in  our  two-plane- wave  model,  the  ac¬ 
tual  Fermi  surface  would  have  distortions  re¬ 
sulting  from  14  Bragg  planes).  We  can  estimate 
the  combined  effects  of  the  necessary  further  re¬ 
daction  in  Fermi-surface  area  as  follows.33  We 
note  that,  with  the  trial  function  given  by  (4),  the 
double  surface  integral  in  (1)  takes  the  form 


Iff 


where  the  rest  of  the  integrand  has  a  somewhat 
weaker  dependence  on  the  location  on  the  Fermi 
surface.  The  ratio  of 


to  a  corresponding  quantity  for  the  free-electron 
case  is  thus  an  approximation  of  the  correction 
factor  required.  Using  inversion  symmetry,  we 
have 

f  (ds^vj  f  {ds2'i,)  |v,  -  v,  I3 
J  (ds?/v?)J  (ds fTrf)  >v“-v°  i3 

f(ds/v)  [dsv  ..  . 

J  (ds°/v°)  J  ds°v°  •  K  ’ 

where  / ds  v  'fds 0  v°  is  just  m/inrpt  while  ( Jds/v)/ 

( f  ds°/v°)  is  commonly  referred  to  as  the  specific- 
heat  effective  mass  »i a/»i.  Combining  this  with 
the  factor  (wop../»i)3,  we  arrive  at  a  total  correction 
factor  /  =  (»iopI/m)(md>/»i).  In  the  two-plane-wave 
model  nzopl  and  can  be  obtained  without  further 
approximation  in  closed  form6: 

where 

a)o=2[S*S/.+  V*(K)],/1-  SK  (19) 

and 

w,  =  2[SyS  SK. 

Further, 

~ =  1  "  \  f  (1  -  x~  1  +  +  «  w3(A)]l/3}l/3) , 

where  X=K'2kf.  These  formulas  then  serve  as  an 
indication  of  the  accuracy  of  our  numerical  proce¬ 
dure  within  the  two-plane-wave  model.  The  re¬ 
sults  are  listed  in  Table  I.  We  note  also  that  the 
scaled  correction  factors,  i.e. ,/(>•,)  /(rj0) ,  com¬ 
pare  very  well.  The  resistivity  of  aluminum  as  a 
function  of  compressed  volume  is  plotted  in  Fig.  3. 
Here  we  display  the  results  using  both  the  correc¬ 
tion  factors  from  the  closed  form  solutions  (/,) 
and  those  from  the  numerical  calculations  (/„). 

The  equation  of  state  of  aluminum  has  been  ob¬ 
tained  by  Friedli  and  Ashcroft.’  They  examined 
most  of  the  common  crystal  structures  and  con¬ 
cluded  that  for  the  pressure  range  they  considered 
(up  to  and  above  300  GPa)  the  fee  structure  is  that 
of  lowest  calculated  energy.  We  can  numerically 
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TABLE  I.  Quantities  mopt/"i  and  evaluated  in  the  two-plane-wave  model.  Both 

closed-form  (analytical)  and  numerical  results  are  listed.  The  quantity/  Is  defined  by 
/-  <>nop,/»i)(m,h/»i).  Wote:  r,6= 2.073.) 


/(rt)//(rl0) 


V/V, 

Analytical 

Numerical 

Analytical 

Numerical 

Analytical 

Numerical 

2.073 

1.00 

1.370 

1.305 

0.985 

0.987 

1.00 

1.00 

2.001 

0.90 

1.644 

1.S61 

0.976 

0.979 

1.19 

1.18 

1.924 

0.60 

2.070 

1.910 

0.964 

0.966 

1.48 

1.43 

1.S41 

0.70 

2.819 

2.552 

0.947 

0.949 

1.98 

1.88 

eliminate  r,  between  this  equation  of  state  and  the 
resistivity  variation  we  obtain  in  this  calculation 
to  arrive  at  a  scaled  resistivity  curve  which  is 
shown  in  Fig.  4.  Unfortunately,  we  have  not  been 
able  to  locate  experimental  data  on  the  pressure 
variation  of  the  resistivity  of  crystalline  alumi¬ 
num.  The  data  of  Bridgman  summarize  the  mea¬ 
surement  of  relative  resistance  rather  than  the 
relative  resistivity.23'3'*  However,  we  can  get  an 
estimate  of  the  experimental  value  of 

[A(p/po)/A(  V/  V()]r/r0*i 

from  the  Bridgman  data  by  using  the  approximate 
equation  for  an  isotropic  cubic  crystal  (roughly 
applicable  to  the  experimental  arrangement  of 


FIG.  3.  Scaled  resistivity  of  aluminum  as  a  function  of 
compressed  volume  at  T  =  300 ‘K  In  the  single-Bragg  - 
planc  two-plane-wave  approximation.  The  curve  labeled 
f,  is  obtained  by  applying  the  Fermi-surface-avea  cor¬ 
rection  factor  in  closed  form:  the/,  curve  is  obtained 
with  the  correction  factor  obtained  numerically  (see 
text). 


Bridgman): 

*(p/Po)  I  »  (20) 

MV/VJ  I  K/,vl  &[V/Vf)  I  K/r0»,  ^ 

With  this  equation  we  get  2.3  for  the  experimental 
value*5  of  A(p/p0)/A(  U/vy  and  this  should  be 
compared  with  the  theoretical  result  of  2.5  ob¬ 
tained  here. 

VI.  DISCUSSION  AND  CONCLUSION 

Our  calculations,  though  they  describe  quite  well 
the  qualitative  trend  of  the  behavior  of  resistivity 
under  pressure  (as  compared  with  experiments  at 
low  pressures),  do  not  yield  particularly  accurate 
numerical  results  for  the  resistivity  itself.  At 
normal  conditions  the  resistivity  p(0)  is  overesti¬ 
mated  by  50'ji,  (the  computed  values  are  4.2  and  4.0 


FIG.  4.  Scaled  resistivity  of  aluminum  as  a  function 
of  pressure.  The  curves  are  obtained  by  numerically 
eliminating  the  volume  between  the  curves  shown  in  Fig. 
3  with  the  equation  of  state  of  aluminum.  The  equation 
of  state  is  obtained  by  the  procedure  described  in  Ref. 

3. 
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r if! cm  with  the  correction  factors/,  and /„  re¬ 
spectively)  and  is  characteristic  of  the  errors  en¬ 
countered  in  the  equivalent  calculation  in  the  liquid 
state.1  The  discrepancy  can  be  attributed  to  the 
approximations  we  have  made  and  to  the  choice  of 
trial  function.  The  Fermi  surface  over  which  we 
perform  our  double  surface  integral  is  at  best  a 
rough  approximation  of  the  real  Fermi  surface 
which  has  a  much  more  complex  geometry.  Our 
approximation  of  a  single  Bragg  plane  allows  for 
more  transitions  to  be  included  as  the  Fermi  sur¬ 
face  thus  generated  has  a  larger  area  than  the  real 
Fermi  surface.  This  additional  surface  is  only 
approximately  accounted  for  in  our  treatment  of 
applying  the  correction  factor  (18).  Another  source 
of  error  can  be  traced  to  the  replacement  of  the 
phonon  frequencies  by  three  branches  (one  longi¬ 
tudinal  and  two  identical  transverse)  and  the  re¬ 
placement  of  the  first  Brillouin  zone  by  a  Debye 
sphere.  This  simplifies  the  problem  but  at  the 
same  time  discounts  all  the  effects  of  phonon  an¬ 
isotropy.  The  errors  introduced  by  using  a  single 
Griineisen  parameter  to  describe  the  phonon  fre¬ 
quency  changes  will  be  averaged  out  to  a  certain 
extent  (at  least  for  small  volume  changes)  as  an 
integral  is  performed  to  arrive  at  the  resistivity. 
Furthermore,  we  note  that  the  two-plane-wave 
approximation  is  not  adequate  for  some  electronic 
levels,  especially  those  near  the  zone  edges  and 
corners.  The  electronic  levels  near  the  symme¬ 
try  point  IF,  for  instance,  require  four  plane  waves 
for  an  adequate  characterization.  A  realistic  de¬ 
scription  of  the  electronic  levels  is  attempted  only 
near  the  Bragg  plane  under  consideration— and 
this  is  only  done  crudely— a  simplification  '  at  in¬ 
evitably  introduces  errors  into  our  results.  We 
note,  however,  that  the  fraction  of  the  Fermi  sur¬ 
face  that  is  close  to  a  Bragg  plane  is  already 
small.  Thus,  we  can  expect  that  the  fraction  that 
lies  near  the  intersection  of  two  or  three  of  these 
Bragg  planes  to  be  smaller  still.  Further  Bragg 
planes  [the  next  class  is  (220)]  are  excluded  from 
our  consideration  as  they  lie  outside  the  Fermi 
sphere.  The  inadequacy  of  the  trial  function  we 
use  necessarily  leads  to  overestimation  of  the  re¬ 
sistivity. 

However,  we  believe  our  results  for  the  sealed 
quantity  have  quantitative  validity  in  the  descrip¬ 
tion  they  give  of  the  changes  in  the  resistivity  as 
pressure  is  applied,  the  previous  objections  not¬ 
withstanding.  The  reason  is  that  when  we  consider 
the  changes  brought  about  by  the  application  of 
pressure  then  the  relative  uncertainties  in  each  of 
the  considerations  should  bo  considerably  re¬ 
duced. 

The  results  of  our  calculation  are  somewhat  in¬ 
teresting  in  that  the  resistivity  of  aluminum  shows 


a  minimum  as  the  metal  is  compressed.  The 
mechanisms  that  cause  this  minimum  are  not  the 
same  as  those  believed  to  account  for  the  resistiv¬ 
ity  minima  in  the  alkali  metals.  In  the  latter  case 
the  admixture  of  d  states  in  the  levels  at  the  Fer¬ 
mi  energy  is  believed  to  be  responsible, **.*•-» 
while  in  the  case  of  aluminum  the  minimum  ap¬ 
pears  to  be  a  manifestation  of  changes  in  the  Fer- 
mic-surface  geometry.  In  our  simple  model  of 
the  electron-ion  pseudopotential  under  compres¬ 
sion  we  see  that  for  all  the  compressed  volumes 
considered  the  reduced  Fourier  coefficients 
£.‘#(lll)/5,  and  vt(200)/Sf  increase  as  the  metal  is 
compressed.  These  increases  cause  the  reduc¬ 
tion  in  Fermi-surface  area  (as  compared  to  the 
free-electron  sphere)  to  increase  from  about  14% 
to  about  38%  when  the  volume  of  the  sample  is  re¬ 
duced  by  30%.  The  values  of  the  reduced  Fourier 
coefficients  and  the  Fermi  surface  area  reduc¬ 
tions  are  listed  in  Table  II. 

It  is  quite  clear  that  a  more  careful  treatment 
of  the  resistivity  variation  requires  that  the  Fer¬ 
mi-surface  distortions  be  fully  taken  into  account 
If  the  Fermi-surface  effects  are  included  (as  de¬ 
scribed  in  Sec.  V),  we  see  that  the  simple  picture 
of  a  one-plane -wave  treatment  needs  to  be  modi¬ 
fied.  In  addition  to  the  (main)  effect  of  reductions 
in  resistivity  stemming  from  increases  in  phonon 
frequencies21  we  have  an  offsetting  effect  from  the 
increase  in  the  distortions  of  the  Fermi  surface 
which  ultimately  reverses  the  trend.  The  net  cor¬ 
rection  factor  that  should  be  applied  for  these  dis¬ 
tortions  should  be 

/=  &»oP|/»l)('«rt/«l) 

instead  of  the  (m0~, An)*  appearing  in  Eq.  (1)  be¬ 
cause  of  a  corresponding  reduction  in  the  region 
of  integration  for  the  double  surface  integral  in 
Eq.  (1).  As  nt* An  remains  near  unity  for  all  the 
values  of  r,  considered,  we  see  that  the  correction 
factor  is  reflectedby  the  increase  (as  a  function  of 
compressi  on)  of  the  optical  mass.  The  two  effects 
together  cause  a  resistivity  minimum  at  !’  F0 
so. 8  which  would  be  absent  in  the  one-plane-wave 
treatment. 

TABLE  II.  Reduced  Fourier  coefficient  r,U ll)/§f r  and 
r,(200)/|cf  and  t-e  reduction  in  Fermi-surface  area  as  a 
function  of  compressed  volume.  The  two-plane-wave 
model  Is  used  to  compute  the  Fermi-surface  (FSi  area 
reductions. 


’V"„> 

I'/F, 

.■ylll)/?«. 

iy  200) 

FS  Reduction 

2.0:3 

1.00 

0.0080 

0.0949 

14  f7 

2.001 

0.90 

0.0322 

0.1144 

22* 

1.924 

0.80 

0.(  585 

0.1350 

29  <7- 

1.84 

0.70 

0.0:70 
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While  the  validity  of  the  simple  one -plane -wave 
treatment  can  be  justified  for  high  temperatures 
and  zero  pressure,  we  have  found  that  in  dealing 
with  the  compressed  metal  a  more  realistic  treat¬ 
ment  taking  into  consideration  the  Fermi-surface 
distortions  is  plainly  required.  As  discussed 
above,  this  differs  from  the  situation  in  the  alkali 
metals  where  the  distortions  in  the  Fermi  surface 
are  second-order  effects.  As  another  example  of 
a  polyvalent  metal  we  may  consider  Pb  whose 
Fourier  coefficients  (or  rather  their  magnitudes) 
are  expected  to  decrease  as  the  metal  is  com¬ 
pressed.30  Therefore,  both  the  Fermi-surface 
distortion  (as  reflected  in  the  optical  mass)  and 
the  increase  in  phonon  frequencies  cause  the  re¬ 
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We  calculate  the  zero-temperature  equation  of  state  of  potassium  with  a  model  Hamil¬ 
tonian  that  includes  core-polarization  effects.  Density  fluctuations  in  the  ion  cores  lead  to 
van  der  Waals  interactions  that  are  dynamically  screened  by  the  valence  electrons.  They 
also  lead  to  screening  of  other  static  interactions,  effects  that  are  incorporated  through  the 
use  of  a  background  dielectric  function  ec(q).  Inclusion  of  core-polarization  effects  yields 
significant  improvement  between  the  theoretical  and  experimental  equations  of  state,  partic¬ 
ularly  at  high  pressures. 


I.  INTRODUCTION 

In  calculations  of  the  thermodynamic  functions 
of  the  simple  metals  it  is  common  to  assume  that 
core-polarization  effects  can  be  neglected,  an  as¬ 
sumption  that  is  normally  justified  on  the  grounds 
that  the  electrons  in  the  corresponding  ions  are 
tightly  bound  and  in  consequence  not  significantly 
polarizable.  It  then  follows  that  the  dispersion 
forces  and  the  effects  of  any  background  dielectric 
shielding  associated  with  the  internal  structure  of 
the  ions  should  be  small,  at  least  in  comparison 
with  the  much  stronger  Coulomb  interactions  and 
electron-gas  properties.  Under  normal  conditions 
these  assumptions  are  usually  valid:  The  screened 
van  der  Waals  interactions,  for  example,  produce 
only  a  slight  softening  of  the  replusive  part  of  the 
pseudopotential-derived  ion-ion  potential.1  At  high 
pressures,  however,  where  interionic  separations  are 
much  reduced,  the  softening  of  such  potentials  can 
have  noticeable  effects,  as  we  shall  see  here.  In  ad¬ 
dition,  the  likely  importance  of  the  background 
screening  can  be  gauged  by  examining  the  quantity 
6f(0)  ~  1  +  4-n,a,  where  a  is  an  ionic  polarizabil¬ 
ity  and  n,  is  the  ionic  number  density.  In  the  alkali 
metal  series  this  does  not  depart  appreciably  from 
unity  (the  range  is  about  1.01  — 1.27);  nevertheless, 
such  departures  can  lead  to  quite  significant  correc¬ 
tions  to  the  various  terms  comprising  the  total  ener¬ 
gy  of  a  metal,2  particularly  at  high  pressure  where 
ft,  has  been  increased. 

In  a  previous  paper3  (referred  to  as  I)  we  derived 
a  model  Hamiltonian  to  treat  the  problem  of  a  sys¬ 
tem  of  interacting  dipoles  and  electrons  for  use  in 
the  calculation  of  thermodynamic  and  structural 
properties  of  simple  metals.  It  was  shown  there 
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that  the  collective  core-core  excitations  give  rise  to 
screened  polarization  waves,  and  the  lowest-order 
dispersion  forces  associated  with  these  waves  are  the 
screened  van  der  Waals  interactions.  It  was  also 
shown  that  in  addition  to  participating  in  physical 
processes  involving  collective  excitation  and  dynam¬ 
ic  screening,  the  polarizable  ions  screen  all  the  static 
interactions.  At  the  microscopic  level  the  system  of 
ions  constitutes  an  inhomogeneous  dielectric  which 
for  long-ranged  interactions  in  situations  of  high 
symmetry  may  be  reasonably  well  approximated  by 
a  dielectric  continuum  with  dielectric  function  ec(q). 
In  this  paper  we  shall  apply  the  principal  results  of 
I  to  be  the  case  of  potassium  which  has  a  rather 
substantial  ionic  polarizability  (see  Table  II  of  Ref. 

1).  It  also  has  a  high  compressibility  which  makes 
it  well  suited  to  the  present  calculation. 

The  paper  is  orgainized  as  follows:  In  Sec.  II,  the 
key  results  for  the  total  energy  of  a  simple  metal 
with  polarizable  ions  are  restated.  In  Sec.  Ill  we  in- 
torduce  the  additional  approximation  needed  to 
evaluate  the  zero-temperature  isotherm.  This  re¬ 
quires  us  to  address  the  form  of  the  static  back¬ 
ground  dielectric  function  ec{q)  for  a  dielectric  con¬ 
tinuum  representing  the  ions,  as  well  as  the  atomic 
polarizabilities  themselves.  The  results  are  discussed 
in  Sec.  IV. 

At  high  compression  the  effects  we  are  discussing 
are  quite  significant,  as  will  be  seen.  It  should  be 
noted  here  that  by  ignoring  core-polarization  effects 
entirely  it  is  still  quite  possible  to  obtain  a  detailed 
quantitative  description  of  the  elastic  properties  of 
the  alkali  metals4,5  using  pseudopotential  methods. 

We  shall  see  below  that  the  explanation  for  this 
minor  paradox  is  simply  the  observation  that  for 
modest  compressions  the  correction  terms  arising 

1636  ©1981  The  American  Physical  Society 


24 


CORE  POLARIZATION  AND  THE  EQUATION  OF  STATE  OF  . . . 


1637 


from  core  polarization  effects  can  be  almost  entirely 
absorbed  into  the  standard  zero  pressure  fitting  pro¬ 
cedure  in  which  the  aggregate  of  all  long-wavelength 
terms  is  adjusted  to  an  experimental  datum.  The 
volume  dependence  of  the  corc-polarization  terms 
is,  however,  sufficiently  different  that  under  more 
extreme  circumstances,  such  as  the  treatment  of 
shock  Hugoniots  and  the  determination  of 
Gruneisen  parameters,  this  procedure  may  require 
reexamination. 


II.  GROUND-STATE  ENERGY 
OF  A  SIMPLE  METAL 


Wc  shall  neglect  the  energy  associated  with  the 
nuclear  degrees  of  freedom5:  The  Hamiltonian  for  a 
simple  metal  of  valence  Z  and  of  volume  fl  whose 
ions  are  polarizable  can  then  be  written1'1 


+  2  j* W  -  q> 


(la) 

(lb) 


,  .  v  2\(q) 


-V] 


(1c) 


As  described  in  I  the  terms  here  have  the  following 
meaning:  Term  (la)  is  the  Hamiltonian  of  an  in¬ 
teracting  electron  gas  in  a  neutralizing  dielectric 

continuum  with  dielectric  constant  ee(q).  The _ 

valence-electron  density  operator  is  p"(  <))=£, e"1  r<; 
the  Coulomb  interaction  is  ve(q)  =  4ire2/q2. 

Term  (lb)  is  a  characteristic  pseudopotential  form 
for  the  interaction  between  the  valence  electrons  and 
the  ions,  but  modified  here  to  incorporate  the  effects 
of  the  dielectric  continuum.  The  pseudopotential 
VjJq)  used  in  (lb)  is  assumed  local;  notejhat  the 
ion  density  operator  is  pVq)  =  Re,qR.  Term 

(lc)  is  the  Coulomb  energy  of  point  charges,  also  in 
the  dielectric  continuum.  The  terms  in  (Id)  are, 
respectively,  the  sum  of  all  q  =  0  terms  and  the 
screened  van  der  Waals  interactions.  As  is  well 
known,  £0  has  an  inverse  volume  dependence6:  If 
rt,  is  the  mean  valence  electron  density,  it  can  be 
written  as 

Ea  —  — — 

477^/3 

where 

r,a0  =  (3/4im,),/J  . 


+  SZE0 + j 


+ 


R  -  R) 


(Id) 


The  screened  fluctuating  dipole  interaction  between 
ions  is1,7 


**(/•) 


d2ufr,iu) 

1 

2 

6u(r,/u) 

l 

dr 2 

dr 

(3) 


where 


fn(q,iu) 


and 


Em- 


the  quantity  etl(q,eo)  being  the  wave-number-  and 
frequency-dependent  dielectric  function  of  the  in¬ 
teracting  electron  gas. 

Perturbation  theory  can  be  applied  in  a  straight¬ 
forward  manner  to  (1).  The  contributions  to  the 
ground-state  energy  are  then:  (a)  the  corrected  ener¬ 
gy  of  the  interacting  electron  gas,  which  in  Ry  per 
electron,  is* 

E'n  =  E,t  +  ,  (4a) 


(4c) 

and  where  Stf{q)  is  the  static  structure  factor  of  the 
interacting  electron  gas;  (b)  the  corrected  Madelung 
energy9 

Em  —  Em  +  Emc  »  (5a) 

where 

aMZxn 

Em  =  -  .  (5b) 

*» 


where 


0.916 

r. 


+  E  t 


Emc  =  T^2  Zvt(q)SMft-  '<*)  -  1  ]  ,  (5c) 

(4b) 

and  S ^  is  the  ensemble  averaged  ionic  structure 
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factor;  (c)  the  corrected  band-structure  energy 


III.  ZERO-TEMPERATURE  ISOTHERM 


E 


m  _  1  NZ 
W  2  fl 


X 


ec(g) 


r(<7)Sij9)  , 


(6) 


where  Xv(q)  is  the  first-order  polarizability  of  the  in¬ 
teracting  electron  gas;  and  (d)  the  aggregate  of  all 
q  =*  0  terms,  as  given  by  (2).  If  these  terms  are  aug¬ 
mented  by  the  screened  fluctuating  dipole  interac¬ 
tions,  the  total  static  lattice  ground-state  energy  can 
be  written 


E(rt)  =  NZ(E'tt  +  Em  +  Eg*  +  E0 ) 

+  7  2  <t>L  (R  -  R'j  •  (7) 

R.R' 


Here  the  band-structure  energy  has  been  calculated 
to  second  order  in  the  electron-ion  pseudopotential. 
For  monovalent  metals  (such  as  potassium)  this 
neglect  of  higher-order  band-structure  energy  terms 
is  satisfactory.10  The  corrections  associated  with 
ec(q)  are  also  calculated  within  a  linear  response 
framework,  again  a  reasonable  approximation  since 
the  effects  we  are  incorporating  here  are,  in  any 
event,  fairly  small. 


J 


If  F  is  the  Helmholtz  free  energy  of  a  system  of 
N  particles  at  temperature  T,  then  the  pressure  is 
p  —  —  (bF  /bfl)TK.  At  T  =  0,  F  =  E,  which  is 
given  by  (7).  The  various  contributions  to  E  require 
the  averaged  core-dielectric  function  f,(g),  the 
frequency-dependent  valence-electron  dielectric 
function  f(g,cu),  the  f-equency-dependent  ionic  po¬ 
larizability  a(<u),  and  the  valence-electron  static 
response  functions. 

Provided  the  background  dielectric  constant 
€c{q  =  0)  is  reasonably  close  to  unity,  as  is  the  case 
here,  we  may  neglect  local-field  corrections  and 
write,  as  earlier 

ec(q)  —  1  +  4n7t,a(0)  (q— *0)  .  (8) 


In  the  opposite  limit  (q  >  2 ir/d,  where  d  is  an  ionic 
diameter)  we  make  a  local-density  approximation, 
taking  the  large  wave-vector  limit  of  the  Lindhard 
result10 


ec(q)  -  1  + 


1 

lTrkFa0 


(9) 


Here  kcF  is  a  characteristic  Fermi  wave  vector  corre¬ 
sponding  to  a  locally  uniform  core-electron  density 
ne.  Let  x  =  (q/7kF)\  then  for  Z  =  1,  a  simple  in¬ 
terpolation  between  (8)  and  (9)  is 


ec(q)  =  1  +  4n7i,a(0) 

1  +  4 

*v_ 

a(0) 

3 

{kFa0)*x* 

»« 

a° 

where  a(0)  =  a(cu  =  0)  and  kf  —  hr2  nv  (u  refers  to  the  valence  electrons).  This  form  for  ec{q)  is  similar  to 
the  wave-vector-dependent  dielectric  functions  used  in  semiconductors  and  other  narrow-gap  insulators."  It  is 
plotted  in  Fig.  1  for  the  choice"  a(0)  =  0.9  A5  and"  ne  =  8 n„. 

To  evaluate  the  correction  to  the  standard  electron-gas  energy,  we  note  that  the  second  term  of  (4)  can  be 
identified  as  one-body  and  two-body  contributions,  that  is  , 


E '  -  —  T  vl(q)— — 


eclq) 


-  1 


et{q) 


-  1 


(11) 

(12) 


For  an  empty-core  pseudopotential14  with  core  radius  re  and  the  interpolation  form  (10),  the  first  of  these  can 
be  evaluated  in  closed  form  and  the  result  is  (in  rydbergs) 

—  (fcffl0)(G/2fl3)[l  +  e-'/3"(cosv/2ai  +  siWlaj)]  ,  (13) 


I - 

with  and 

G  =  (nc/nv)(3kFa0)~'  ,  s  =  lkFrt  . 

a*  =  G[1  +  4/t1a(0)]/4s-nia(0)  ,  To  determine  (12)  it  is  sufficiently  accurate  to  take 
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FIG.  1.  Interpolated  form  for  ec(g)  for  potassium. 
The  dotted  curves  show  the  high-  and  low-?  limits  of 
ec(q).  The  solid  curve  is  from  Eq.  (10).  (a)  r,  =  4.86 
(ft/n0  =  1.0);  (b)  r,  =  4.10  (fl/n0  =  0.6). 

the  Hartree-Fock  result15  for  the  electron-gas  struc¬ 
ture  factor  Sttiq),  i.e., 

S'Jq)  =  yJC  -  y*3  lx  <  1)  , 


Sttlq)~  1  (X  <  1)  , 


where 

x  =  q/lkF  . 

Since  ec(x)  is  already  close  to  unity,  it  is  reasonable 
to  ignore  further  corrections  to  the  correlation  ener¬ 
gy  arising  form  the  modification  of  the  effective 
electron-electron  interaction  because  of  core- 
polarization  effects. 

We  turn  now  to  the  energy  of  the  screened  van 
der  Waals  interactions:  Its  determination  requires 
the  frequency-dependent  atomic  and  electron-gas 
polarizabilities.  Again,  since  this  contribution  is  not 
a  major  one,  we  may  make  reasonable  approxima¬ 
tions  in  both  quantities.  For  the  atomic  polarizabili¬ 
ty  we  take  the  usual  Lorentz  form 

ala)  =  a(0)[l  —  (m/wo)2]"1  ,  (15) 

where  coQ  is  a  characteristic  frequency  (estimated  to 
be  the  equivalent  of  47  eV  for  potassium16).  The 
wave-vector-  and  frequency-dependent  dielectric 
function  for  the  valence  electrons  is  taken  to  have 
the  simple  form17 

clq,(i))  =  1  +  kj^/lq*  —  kjfa^/o)p) ,  (16) 

where  k^  is  the  Thomas-Fermi  wave  vector.  No¬ 
tice  that  this  leads  to  the  familiar  q  -*0  (for  tu  =  0) 
and  fa  — t-0  (for  q  =  0)  limits.  The  point  is  that  the 
use  of  (151  and  (16)  in  conjunction  with  (3)  gives 

NZE^W  =  y  £  *£(R)  , 

R*  0 

where 


rLiR)=- 


e1  3/wq  1 

<*o  2 

ao 

2a<j  4(eJ/2a0)  [1  +  lo)p/a0)3] 

a° 

R 

with 


fieOp 


I 


(e 2/2<j0)  2 ir 


f  dy 


aHiu ) 


a  o 


/i(R.i«) 


(17) 


FtlR,iu )  =  [1  +  { u  /(Op  )2]-  'cosh4Qree  " 1QR 


°o 

6 

ao 

X 

6 

R 

+  12  Qa0 

R 

10«2a0)J 


£o 

R 


+  4(Qa0)J 


£o 

R 


+  «2a<>>* 


fo 

R 


2 


Here  the  quantity  Q  is  defined  by 

Q 2  =  k  2f{!  +  (u/ap)1]  . 

The  sum  in  (17)  can  be  evaluated  with  relatively  lit¬ 
tle  computational  effort.  Finally,  we  determine  the 
parameter  a  in  Eq  by  the  zero-pressure  condition6 


which  leads  to 


a 


9 


r,0(0.916  +  Z2/3cr.v)  +  0.03  lr* 


-  f.0 


Ml 

dr, 


—  4.42  r,0 

f.o 


(18) 
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where  E'  refers  to  the  sum  of  the  energy  contribu¬ 
tions  other  than  those  having  their  origin  in 
electron-gas  and  electrostatic  terms.  In  the  case  of 
potassium  the  zero-pressure  electron  spacing  param¬ 
eter18  is  r, o  =  4.860. 

The  total  energy  of  bcc  potassium  can  now  be  ob¬ 
tained  by  adding  the  contributions  given  by  Eqs. 

(4)  — (6)  and  (17)  with  a  chosen  to  give  the  correct 
zero-pressure  density,  and  subsequently  held  fixed. 
The  various  contributions  to  the  energy  and  pressure 
are  listed  in  Table  I  for  a  sequence  of  volume 
compressions.  The  T  =  0  equation  of  state  is 
shown  in  Fig.  2  and  compared  there  with  the  low- 
temperature  experimental  results  of  Monfort  and 
Swenson.1’  Above  20  kbar  the  comparison  is  made 
to  an  extrapolation  achieved  using19 

Pr  =  0(n>  =  37.0[(lVn)3  85  -  1  J/3.85  (19) 

in  units  of  kbar.  It  is  clear  from  Fig.  2  that  the 
agreement  between  theory  and  experiment  is  satis¬ 
factory. 

The  importance  of  core-polarization  contributions 
can  be  judged  by  comparing  these  results  with  core 
effects  neglected  altogether.  In  practice  this  is 
(unintentionally)  carried  out  by  the  fitting  procedure 
in  which  all  q  =  0  terms  are  assumed  incorporated 
in  the  aggregate  E0.  Since  the  actual  zero  pressure 
density  is  determined  by  all  contributions  to  the  en¬ 
ergy,  the  explicit  neglect  of  core  polarization  is 
equivalent  to  assuming  their  implicit  incorporation 
in  E0  (and  assigning  them  thereby  an  inverse 


volume  dependence).  Figure  3  and  Table  II  sum¬ 
marize  the  results  of  such  a  procedure.  A  direct 
comparison  of  these  two  sets  of  results  show  that 
though  there  is  a  qualitative  similarity  (and  quantita¬ 
tive  to  within  30%)  the  effects  of  core  polarization 
are  not  small  in  potassium.  The  comparison  also 
shows  that  to  the  extent  that  the  difference  in  ener¬ 
gies  computed  by  the  two  different  methods  can  be 
assumed  inversely  volume  dependent,  much  of  the 
difference  can  be  absorbed  in  Eq. 

IV.  DISCUSSION 

As  might  have  been  anticipated  from  the  fact  that 
the  potassium  ion  is  quite  polarizable,  we  find  that 
core-polarization  effects  in  postassium  metal  contri¬ 
bute  noticeably  to  its  thermodynamic  functions,  par¬ 
ticularly  at  high  pressures.  We  now  examine  the 
degree  to  which  this  conclusion  may  depend  on  the 
approximations  and  numerical  procedures  adopted. 
First,  the  interpolation  form  used  for  ee(q)  is  plausi¬ 
ble  on  physical  grounds,  is  convenient  for  numerical 
reasons,  but  is  nevertheless  still  approximate.  On 
the  other  hand,  much  of  the  contribution  from  these 
terms  originates  with  the  small-9  contribution  where 
the  background  dielectric  constant  is  best  known. 

In  the  same  long-wavelength  limit  the  dipole  ap¬ 
proximation  used  throughout  is  expected  to  be  valid 
and  so  is  the  standard  procedure  for  evaluating 
local-field  effects.20  When  wavelengths  become 
comparable  to  the  spatial  extent  of  the  core-electron 


TABLE  I.  Computed  energy  and  pressure  for  potassium  as  a  function  of  volume  with  core-polarization  effects  taken 


into  account. 

The  various 

quantities  are 

defined  in  the 

text.  fl0  is 

the  zero-pressure  volume. 

n/ft0 

Em 

Energy  (Ry/eleetron) 

Eo  Emc 

E»  En 

E  wjw 

em 

1.0 

-0.1610 

-0.3687 

0.2793 

-0.1003 

-0.0149  -0.0557 

-0.0012 

-0.4225 

0.9 

-0.1620 

—  0.3819 

0.3104 

-0.1085 

-0.0171  -0.0609 

-0.0015 

-0.4215 

0.8 

-0.1628 

-0.3972 

0.3492 

-0.1184 

-0.0199  -0.0673 

-0.0019 

-0.4183 

0.7 

-0.1633 

-0.4152 

0.3991 

-0.1302 

-0.0237  -0.0751 

-0.0024 

-0.4108 

06 

—  0.1633 

-0.4371 

0.4656 

-1448 

-0.0287  -0.0852 

-0.0032 

-0.3967 

n/n0  P(E 

r|  +  E»  +  Eq) 

P(Evdw) 

Pressure  (kbar) 

PIEmc)  P(Ek) 

P(E*) 

PM 

P  1 

1.0 

44.6 

-0.7 

-23.3 

-6  0 

-14.5 

0.0 

0.0 

0.9 

59.3 

-1.0 

-27.5 

-7.6 

-17.5 

5.7 

4.8 

0.8 

80.7 

-1.3 

-32.8 

-9.9 

-21.5 

15.3 

13.1 

0.7 

113.0 

-2.0 

-40.0 

-  13.1 

-27.0 

30.9 

(28.3) 

0.6 

164.5 

-3.0 

-49.8 

— 18.1 

-35.2 

58.5 

(59.1) 

'Reference  19.  Numbers  in  parentheses  are  obtained  by  extrapolation. 
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FIG.  2.  Comparison  of  computed  equation  of  state  of 
potassium  (core-polarization  effects  included)  with  the  ex¬ 
perimentally  determined  function.  The  dotted  portion  of 
the  experimental  curve  is  obtained  by  extrapolation  of 
Eq.  (32)  (Ref.  19). 

charge  distribution,  however,  the  description  in 
terms  of  translationally  invariant  core-electron 
response  functions  is  not  strictly  justified.  For  such 
wavelengths  it  is  the  case  that  the  dielectric  function 
approaches  unity  quite  closely  so  that  though  the  in¬ 
terpolation  form  we  use  [Eq.  (10)]  lacks  the  expect¬ 
ed  symmetry,  the  error  thereby  introduced  is  cer¬ 
tainly  small. 

It  is  also  worth  noting  that  the  parameters  a(0)  and 


Mo  can  be  expected  to  be  weak  functions  of  densi¬ 
ty.21  Furthermore,  it  is  well  known  that  a  nonlocal 
pseudopotential  is  required  for  a  correct  decription 
of  band  structure  of  potassium  (there  are  empty  d- 
symmetry  bands  lying  above  the  Fermi  energy22-35). 
Of  the  terms  most  affected  by  such  considerations, 
the  band-structure  energy  is  foremost  Nevertheless, 
once  such  nonlocal  contributions  are  averaged,  as  in 
the  construction  of  the  total  energy,  the  resulting 
uncertainties  in  the  equation  of  state  are  not  great 
The  reason  is  that  the  contributions  to  the  equation 
of  state  from  EB s  are  —  6  kbar  at  fl/flo  =  1,  and 
-  18  kbar  at  fl/flo  =  —  0.6.  If  EBS  were  complete¬ 
ly  omitted,  then  at  p  =  0,  E0  would  acquire  (by  vir¬ 
tue  of  the  fitting  procedure)  an  additional  term  to 
cancel  the  —  6  kbar.  Since  E0  is  assigned  an  inverse 
volume  dependence,  it  would  yield  — 17  kbar  at 
fl/n0  —  0.6,  thus  leading  to  a  total  error  of  only  1 
kbar.  At  very  high  pressures,  the  d  bands  which  are 
the  major  source  of  this  nonlocality,  might  actually 
intersect  the  Fermi  surface,  a  phenomenon  ob¬ 
served24  in  Cs  at  around  42.5  kbar.  The  nearly- 
free-electron  approach,  on  which  our  calculations 
have  been  based,  will  then  be  invalidated.  We  note 
that  in  the  experimentally  determined  equation  of 
State  (at  room  temperature)  no  such  transition  of 
this  nature  has  so  far  been  observed  for  pressure  up 
to  50  kbar.19-25  All  of  these  effects  may  need  ulti- 


6 


0 


TABLE  II.  Computed  energy  and  pressure  for  potassium  as  a  function  of  volume  (core  po¬ 
larization  effects  neglected). 


n/n0 

Energy  (Ry/electron) 

Ev  +  Em  +  Eo  ^bs 

Etat 

1.0 

-0.3702 

-0.0166 

-0.3868 

0.9 

-0.3667 

-0.0195 

-0.3862 

0.8 

-0.3607 

-0.0233 

-0.3840 

0.7 

-0.3507 

-0.0283 

-0.3790 

0.6 

-0.3346 

-0.0348 

-0.3694 

Pressure 

(kbar) 

fl/flo 

PiEn  •+■  Em  +  Eo) 

F(£k  1 

P  * 

•  cipt 

1.0 

7.9 

-7.9 

0.0 

0.0 

0.9 

14.1 

-10.2 

3.9 

4.8 

0.8 

23.4 

-13.3 

10.2 

13.1 

0.7 

38.2 

-17.4 

20.8 

(28.3) 

0.6 

62.7 

-22.8 

39.9 

(59.1) 

•Reference  19.  Numbers  in  parentheses  are  obtained  by  extrapolation. 
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FIG.  3.  The  computed  equation  of  state  shown  in  this 
figure  neglects  core  effects. 
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